Mon. Not. R. Astron. Soc. 000, ITH201 (20071 Printed 8 July 2008 



(MN JATfcjX style file v2.2) 



Resurrecting the Red from the Dead: Optical Properties of 
BCGs in X-ray Luminous Clusters* 



oo 
O 

o 

(N 



Chris Bildfellf, Henk Hoekstra|, Arif Babul, Andisheh Mahdavi 

Department of Physics & Astronomy, University of Victoria, Victoria, BC V8P 1A1, Canada 



Submitted year month day 



3 

OO 

S 1 

6 

CO 

> 
(N 

(N 

o 

oo 
o 



ABSTRACT 

We present measurements of surface brightness and colour profiles for the brightest 
cluster galaxies (BCGs) in a sample of 48 X-ray luminous galaxy clusters. These data 
were obtained as part of the Canadian Cluster Comparison Project (CCCP). The 
Kormendy relation of our BCGs is steeper than that of the local ellipticals, suggesting 
differences in the assembly history of these types of systems. We also find that while 
most BCGs show monotonic colour gradients consistent with a decrease in metallicity 
with radius, 25% of the BCGs show colour profiles that turn bluer towards the centre 
(blue-cores) . We interpret this bluing trend as evidence for recent star formation. The 
excess blue light leads to a typical offset from the red sequence of 0.5 to 1.0 mag in 
(g'-r 1 ), thus affecting optical cluster studies that may reject the BCG based on colour. 
All of the blue-core BCGs are located within ~ lOkpc of the peak in the cluster X-ray 
emission. Furthermore, virtually all of the BCGs with recent star formation are in 
clusters that lie above the L x — T x relation. Based on photometry alone, these findings 
suggest that central star formation is a ubiquitous feature of BCGs in dynamically 
relaxed cool-core clusters. This implies that while AGNs and other heating mechanisms 
are effective at tempering cooling, they do not full compensate for the energy lost via 
radiation. 

Key words: galaxies: clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: evolution - cooling flows - X-rays: galaxies: clusters 



1 INTRODUCTION 

It is well known that the hot intra-cluster medium (ICM) 
in galaxy clusters radiates enormous amounts of energy at 
X-ray wavelengths (L x ~ 10 44 — 10 46 ergs -1 ). Taken at face 
value, this loss of energy from the ICM implies that the gas 
should be cooling, losing pressure support and flowing in- 
ward. Early cooling flow models predicted cooling and star 
formation rates on the order of 300 < M < 1000 M© yr" 1 
(Fabian 1994). The lack of widespread detection of iron lines 
expected from cluster gas cooling below 1-2 keV in XMM- 
Newton observations of cool core clusters (Peterson et al. 
2003) contradicted this model and created the need to incor- 
porate non-gravitational heating mechanisms to counterbal- 
ance the radiative energy loss. While these observations have 
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effectively ruled out accretion flows large enough to form an 
L* galaxy in 100 Myr, they do not preclude that star forma- 
tion is ongoing in cool core clusters at a much reduced rate. 
The level of star formation still allowed is uncertain and it 
is important to investigate this through means independent 
of the X-ray-determined mass flow rate. Of particular in- 
terest is the distinction between models in which heating 
completely balances cooling and those in which heating off- 
sets most but not all of the radiative losses. If significant star 
formation is occurring as a result of a heating imbalance it 
should be reflected in the optical properties of the bright- 
est cluster galaxy (BCG), which is often found deep within 
the cluster potential. A good way to test this is to look for 
links between the BCG optical properties with the X-ray 
properties of their host clusters. In fact, several case stud- 
ies over the last two decades have found evidence for recent 
star formation in the central galaxies of cool core clusters 
(Johnstone et al. 1987, McNamara & O'connell 1989, 1992, 
Allen et al. 1992, Donahue & Voit 1997, Cardiel et al. 1998, 
Crawford et al. 1999, Edge 2001, Edge et al. 2002, Goto 
2005, McNamara et al 2006, Wilman et al. 2006, Edwards 
et al. 2007). Furthermore, the strength of the recent star 
formation indicators is at least weakly correlated with their 
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Table 1. CCCP Targets Observed with MegaCam 



Name z Q-bcg &BCG r l -^off Core R, BCG 

(J2000) (J2000) kpc mag/(") 2 kpc colour kpc Code 



3C295 


0. 


46 


14 


11 


20. 


57 


+52 


12 


09. 


.93 


31 


.911 


.89 
.79 


23. 


283j 


-.055 
- 051 


12.0 


blue 


15 


± 


3 


1 


Abell 115 


0. 


20 


00 


56 


00. 


26 


+26 


20 


32 


.65 


11 


.09^ 


.14 

.12 


23. 


937j 


+)05 
- 005 


143 


red 


_ 






1 


Abell 115N 


0. 


20 


00 


55 


50. 


62 


+26 


24 


37. 


.55 


11 


.941 


.or> 

.07 


22. 


301^ 


-.'005 
006 


10.1 


blue 


11 


± 


1 


1 


Abell 223a 


0. 


21 


01 


38 


02. 


29 


-12 


45 


19. 


.92 


12 


.041 


.4!) 
.42 


23. 


852j 


-.'(122 
- 019 


_ 


red 


_ 






1 


Abell 223b 


0. 


21 


01 


37 


55. 


99 


-12 


49 


10. 


.09 


76 


.661 


.56 
.51 


25. 


267j 


-!oio 

- 010 


8.20 


red 


_ 






1 


Abell 520 


0. 


20 


04 


54 


14. 


06 


+02 


57 


10 


.63 


40 


.36+ 


.37 
.42 


24. 


156^ 


-:oi6 

019 


340 


red 


_ 






2 


Abell 521 


0. 


25 


04 


54 


06. 


89 


-10 


13 


24. 


.61 


37 


.111 


.47 
.51 


23. 


614^ 


-:o2i 

- 025 


32.8 


red 


_ 






1 


Abell 586 


0. 


17 


07 


32 


20. 


31 


+31 


38 


01 


.06 


57 


.961 


.47 

.54 


24. 


038j 


-.013 
- 013 


10.7 


red 


_ 






1 


Abell 611 


0. 


29 


08 


00 


56. 


83 


+36 


03 


23. 


.79 


34 


,40l 


.33 
.28 


23. 


068^ 


-!017 
- 015 


4.0 


red 


_ 






1 


Abell 697 


0. 


.28 


08 


42 


57. 


.58 


+36 


21 


59. 


.38 


70 


.541 


.54 
.47 


24. 


396j 


-!012 
- 010 


20.3 


red 


_ 






1 


Abell 851 


0. 


41 


09 


42 


57. 


.47 


+46 


58 


19. 


.94 


17 


.271 


.20 
.98 


22. 


306j 


-.046 
- 171 


278 


blue 


10 


± 


1 


3 


Abell 959 


0. 


29 


10 


17 


36. 


03 


+59 


34 


01 


.81 


10 


fi7 +0.05 
•"'-.09 


22. 


269j 


-:oi2 

- 021 


36.5 


red 


_ 






2 


Abell 1234 


0. 


166 


11 


22 


29. 


95 


+21 


24 


22. 


.00 


21 


.471 


.24 
.22 


23. 


019^ 


-.'022 

- 020 


_ 


red 


_ 






2 


Abell 1246 


0. 


19 


11 


23 


58. 


83 


+21 


28 


49. 


.84 


55 


.00+ 


.37 
.33 


24. 


490j 


-!012 
- 010 


48.4 


red 


_ 






2 


Abell 1758b 


0. 


28 


13 


32 


38. 


.43 


+50 


33 


36. 


.01 


51 


.531 


.33 
.42 


24. 


227j 


+)11 

014 


25.2 


red 


_ 






1 


Abell 1835 


0. 


25 


14 


01 


02. 


10 


+02 


52 


42. 


.69 


66 


.39+ 


.50 
.51 


24. 


125j 


-!016 
- 015 


6. 1 


blue 


19 


± 


2 


1 


Abell 1914 


0. 


171 


14 


26 


03. 


.90 


+37 


49 


53. 


.57 


19 


.161 


.23 
.28 


24. 


397j 


+)08 
-.009 


85.6 


red 


- 






2 


Abell 1942 


0. 


22 


14 


38 


21. 


88 


+03 


40 


13. 


.34 


43 


.691 


.22 
.22 


23. 


549j 


-.008 
-.009 


4.30 


red 








1 


Abell 2104 


0. 


16 


15 


40 


07. 


94 


-03 


18 


16 


.25 


56 


.481 


.37 
.37 


24. 


244^ 


-.012 
-.013 


7.66 


red 








1 


Abell 2111 


0. 


23 


15 


39 


40. 


52 


+34 


25 


27. 


.46 


26 


.421 


.35 
.31 


23. 


206j 


-.025 
-.02.3 


129 


red 








2 


Abell 2163 


0. 


20 


16 


15 


48. 


99 


-06 


08 


41. 


.21 


26 


.331 


.56 
.47 


23. 


148j 


-.0.36 
-.031 


160 


red 








1 


Abell 2204 


0. 


15 


16 


32 


16. 


98 


+05 


34 


32. 


.84 


130.13; 


f .49 
-.33 


25. 


006j 


-.006 
-.004 


1.15 


blue 


11 


± 


3 


1 


Abell 2259 


0. 


.16 


17 


20 


09. 


66 


+27 


40 


08. 


.29 


44 


.94+ 


.51 
.42 


23. 


599j 


-.021 
-.017 


78.2 


red 








1 


Abell 2261 


0. 


22 


17 


22 


27. 


23 


+32 


07 


57. 


.72 


36 


.101 


.14 

.14 


22. 


897j 


-.007 
-.008 


0.39 


red 








1 


Abell 2537 


0. 


30 


23 


08 


22. 


22 


-02 


11 


31. 


.74 


46 


.371 


.42 

.47 


23. 


605j 


-.016 
-.018 


17.0 


red 








1 


CL0910+41 


0. 


44 


09 


13 


45. 


52 


+40 


56 


28. 


.54 


32 


.641 


.49 
.42 


23. 


358j 


-.0211 
-.025 


3.87 


blue 


40 


± 


8 


1 


CL1938+54 


0. 


26 


19 


38 


18. 


12 


+54 


09 


10 


.19 


148.42; 


fl.7 
-1.4 


25. 


298j 


-.019 
-.016 




red 








1 


MS0440+02 


0. 


.19 


01 


43 


09. 


92 


+02 


10 


19. 


.33 


102.00; 


f .33 
-.28 


24. 


458j 


-.005 
-.004 


0.87 


red 








1 


MS0451-03 


0. 


54 


04 


54 


10. 


.84 


-03 


00 


51 


.39 


15 


.55+ 


.29 
.58 


24. 


418j 


-.015 
-.024 


28.4 


blue 


20 


± 


2 


1 


MS1008-12 


0. 


30 


10 


10 


32. 


34 


-12 


39 


52 


.77 


33 


.10+ 


1.5 
1.4 


23. 


214j 


-.073 
-.072 


9.88 


red 








1 


RXJ1347.5-1145 


0. 


45 


13 


47 


30. 


65 


-11 


45 


09. 


.09 


29 


.611 


.16 
.16 


23. 


154j 


-.012 
-.014 


1.64 


blue 


20 


± 


2 


1 



BCGs (by cluster) that were observed with MegaCam. Spectroscopically determined redshifts (z) of the host clusters are shown. The 
coordinates of the selected BCG are given as (ageg, &bcg)\ see section|A]in the appendix for details. The r' best fit surface brightness 
profile parameters r T e and fi^ are given. The value of fi^ has been corrected for galactic extinction and cosmological dimming. The 
projected offset of the BCG with respect to the peak in host cluster X-ray emission is given as R g. Core colour is given as "red" or 
" blue" based on the shape of the inner colour profile (see section |4j . If a blue core is present the size of the blue region is given as R+ . 
Selection code refers to the degree of ambiguity assessed to the BCG selection; [1] obvious (cD), [2] non-trivial (X-ray necessary), [3] 
highly ambiguous. 



cooling-flow derived mass drop out rates (Johnstone et al. 
1987, Heckman et al. 1989, McNamara 1997). 

However, despite the compelling evidence for recent star 
formation in the BCGs of cool core clusters, the source of 
the cold gas required to fuel this star formation is yet uncer- 
tain. Two competing explanations are that it is a product 
of radiative cooling of the ICM (ie. cooling flow) or that the 
cold gas is deposited during a merger event (c.f. Poole et al. 
2006). The uncertainty surrounding the origin of the cold gas 
persists mainly because previous work focused on samples 
of cool core clusters. However, to resolve this issue requires 
a comparison of BCGs across a wide variety of host cluster 
X-ray morphologies (ie. cool cores and non-cool cores). We 
therefore undertake a comprehensive study of BCGs using a 
cluster sample that spans the full observed range of scatter 
in the cluster X-ray scaling relations (ie. mass and X-ray 
luminosity vs. X-ray temperature). 



We present here an analysis of the surface brightness 
and colour profiles for 48 BCGs. We compare these optical 
data to the global X-ray luminosity (L x ) and temperature 
(T x ) of their host clusters. We aim to discover how the recent 
star formation in BCGs is linked to the balance between 
heating and cooling clusters. This information is crucial to 
our understanding of how baryonic feedback influences the 
evolution of the ICM. 

The structure of the paper is as follows. We describe 
the data in section [2] Section [3] contains a discussion of the 
BCG surface brightness profiles. In section|3]we discuss BCG 
colour profiles. Section [5] gives a comparison of BCG optical 
properties to global cluster X-ray properties. The conclu- 
sions are contained in section [5] The cosmological parame- 
ters used throughout this work are Ho = 70 km s _1 Mpc -1 , 
Q m = 0.3 and fi A = 0.7. 
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Table 2. CCCP Targets Observed with CFH12K 



Name 


z 


"SCO 


&BCG 




... -fi' 

' e 






Rofi 


Core 


R» 






BCG 






(J2000.0) 


(J2000.I 


3) 


kpc 


mag/(") 


2 


kpc 


colour 


kpc 




Code 


Abell 68 


0.255 


00 


37 


06.85 


+09 


09 


24.51 


85 76+'™ 


24.780+ 


.015 
.Dili 


15.4 


blue 


16 


± 


3 


1 


AKwll 90Q 

/\oeii zuy 


u. zuo 


01 


31 


52.54 


-13 


36 


ah nn 

4U.UU 


on 7oe + -5C> 

Oi). 1 OO 


9A 3QA~f 


.011 
.011!) 


1 fi A 
ID. 4 


red 


- 






i 
i 


AKwll 9fi7 
/\Deil ZO / 


n 93 

U. ZO 


01 


52 


41.95 


+01 


00 


9^ 8Q 

zo.oy 


/o.oy_ yQ 


93 QQQ + 

zo. oyy _ 


.017 
.010 


77 ^ 
/ 1 .O 


red 


- 






1 
i 


/\dcii o ( ua 


U.o ( 


02 


.39 


52.70 


-01 


34 


1 Q £9 
lo.OZ 


no. y_ 2 3 


9/1 fi33 + 
Z4.0oo_ 


.012.1 
.027 


1 03 

lyo 


blue 


20 


± 


z 


Z 


AKoll 370k 
/\oeii o ( uo 


n 37^ 

U.O / 


02 


39 


53.10 


-01 


34 


^ 71 
00. / 1 


90 AQ+-62 
2».48_ 62 


99 SI 9 + 
ZZ.OIZ 


.042 
.042 


93 1 
Zo. 1 


red 


- 






9 
Z 


Abell 383 


0.187 


02 


48 


03.38 


-03 


31 


44.93 


26 95+'?? 


22.502+ 


.011 
.010 


0.71 


blue 


28 


± 


5 


1 


Abell 963 


0.206 


10 


17 


03.63 


+39 


02 


49.67 


60 00 + '?5 

""'""-.35 


23.679+ 


.010 
.010 


6.33 


red 


- 






1 


AK^ll 1 fiSQ 


n 1 839 

U. lOOZ 


13 


11 


29.52 


-01 


20 


97 Sfi 
Z / .OO 


Afi RA+ 3 " 
40.04_ 27 


93 Q/1Q + 


.014 
.011 


^ ni 

O.Ul 


red 


- 






1 
1 


Abell 1763 


0.223 


13 


35 


20.12 


+41 


00 


04.30 


24 47+ i? 


22.548+ 


.014 
.015 


7.40 


rod 


- 






1 


Abell 2218 


0.1756 


16 


35 


49.26 


+66 


12 


44.56 


46 34+' 39 


23.602+ 


.014 
.017 


59.5 


rod 








1 


A Koll 991 Q 

/\oeii zziy 


n 99^fi 

U.ZZOO 


16 


40 


19.85 


+46 


42 


,41 30 

4:1. OU 


fi1 £^ + '51 


93 Q/M + 

zo.y4± _ 


.014 
.014 


8 97 
O.Z / 


red 








1 
1 


Abell 2390 


0.2280 


21 


53 


36.84 


+17 


41 


44.10 


18 42+'™ 


22.540+ 


.020 
.020 


3.63 


blue 


25 


± 


5 


1 


CL0024+16 


0.39 


00 


26 


35.68 


+17 


09 


43.48 


60.1±J;| 


23.8361 


.042 
.044 


24.5 


red 








3 


MS0016+16 


0.5465 


00 


18 


33.56 


+16 


26 


16.30 


42.38ti2 


23.797+ 


.044 
.048 


41.2 










3 


MS0906+11 


0.1704 


09 


09 


12.76 


+10 


58 


29.12 


55.35 


23.8091 


.012 
.011 


2.55 


red 








2 


MS1224+20 


0.3255 


12 


27 


13.48 


+19 


50 


56.10 


36.57+'™ 


23.3791 


.0:50 
.OO* 




red 








1 


MS1231+15a 


0.2353 


12 


33 


54.72 


+15 


26 


19.96 


51.931^ 


24.3591 


.015 
.019 




red 








3 


MS1231+15b 


0.2353 


12 


33 


55.35 


+15 


25 


59.11 


50.97^42 


24.1001 


.012 
.014 




red 








3 


MS1358+62 


0.3290 


13 


59 


50.62 


+62 


31 


05.09 


48.32+;™ 


24.0711 


.020 
,o:i2 


4.41 


red 








1 


MS1455+22 


0.2568 


11 


57 


15.12 


+22 


20 


34.48 


24.80t;« 


22.2201 


.000 
.038 


3.29 


blue 


29 


± 


3 


1 


MS1512+36 


0.3727 


15 


11 


22.51 


+36 


36 


21.30 


57.56±;*| 


24.111+. 


.024 
.024 


5.74 


blue 


15 


± 


2 


1 


MS1621+26a 


0.4275 


16 


23 


35.17 


+26 


34 


28.40 


116. l+3 7 


26.0041 


.000 
.042 


74.7 










2 


MS1621+26b 


0.4275 


16 


23 


35.45 


+26 


34 


14.50 


46.51^ 


24.1931 


.047 


41.2 


blue 


21 


± 


5 


2 



BCGs (by cluster) that were observed with CFH12K. Spectroscopically determined redshifts (z) of the host clusters are shown. The 
coordinates of the selected BCG are given as (agee Sbcg)] see section|A]in the appendix for details. The r' best fit surface brightness 
profile parameters j-J* and ^ are given. The value of ^ has been corrected for galactic extinction and cosmological dimming. The 
projected offset of the BCG with respect to the peak in host cluster X-ray emission is given as i? ff. Core colour is given as "red" or 
"blue" based on the shape of the inner colour profile (see section [4} . If a blue core is present the size of the blue region is given as R, . 
Selection code refers to the degree of ambiguity assessed to the BCG selection; [1] obvious (cD), [2] non-trivial (X-ray necessary), [3] 
highly ambiguous. 



2 DATA 

The data studied here are part of the Canadian Cluster 
Comparison Project (CCCPfl The CCCP is a study of an 
X-ray selected sample of 50 clusters from the ASCA catalog 
of Horner (2001) with redshifts in the range 0.15 ^ z ^ 0.55 
and most having X-ray temperatures T x > 5 keV. The 
CCCP is selected such that it spans the full range of ob- 
served scatter in the cluster X-ray and Sunyaev-Zel'dovich 
(SZ) scaling relations. This was done to create maximal 
leverage for the investigation of baryonic feedback effects. 
The combination of X-ray and deep, multi-filter optical 
data makes the CCCP an excellent sample for studying 
the links between BCGs and their host cluster environ- 
ments. The deep optical imaging data was obtained using 
the Canada- France-Hawaii Telescope (CFHT) . The full sam- 
ple comprises a set of 30 clusters observed in g and r' using 
the MegaCam detector and 20 clusters observed in the B 
and R filters using the CFH12K detector which were found 

1 Based on observations obtained with MegaPrime/MegaCam, 
a joint project of CFHT and CEA/DAPNIA, at the Canada- 
France-Hawaii Telescope (CFHT) which is operated by the Na- 
tional Research Council (NRC) of Canada, the Institute National 
des Sciences de l'Univers of the Centre National de la Recherche 
Scientifique of France, and the University of Hawaii. 



in the CFHT data archive at the Canadian Astronomical 
Data Centre (CADC) (see Hoekstra (2007) for more details). 
MegaCam and CFH12K have 0.186"/pix and 0.206"/pix im- 
age scales respectively and with the median seeing of 0.7" 
at the CFHT prime focus these cameras image the central 
r < 2>5hi kpc regions of the target BCGs (~ 10" at a typi- 
cal redshift of the clusters in our sample) over several resolu- 
tion elements (see figure ICl)) . All of the data are taken with 
between 0.5" and 1.0" seeing. The FWHM of the PSFs for 
the individual images are shown on the surface brightness 
and colour profiles in appendix [C] 

Table [T] lists all targets observed with MegaCam and 
Table lists the targets observed with CFH12k. For the 
clusters imaged with MegaCam, the ones with z < 0.4 were 
observed for 4 x 450 seconds in g and 8 x 600 seconds in 
r' while clusters with z > 0.4 were observed for 4 x 600 
seconds in g' and 8 x 600 seconds in r' . The clusters that 
were imaged with CFH12K were taken with a variety of 
exposure times, ranging from 3000 to 15300 seconds in B 
and from 3000 to 13800 seconds in R. The exact exposure 
times for the individual CFH12K clusters can be found in 
Hoekstra (2007). 

The initial data reduction (bias subtraction and flat- 
fielding) are done on site at CFHT with the Elixir auto- 
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mated instrument and software pipelinfl The pipeline also 
provides photometric zero-points based on observations of 
photometric standard stars during the observing run. The 
images of each cluster are aligned and median-stacked using 
the SWarp software package. 

The BCG selection is carried out using the combined 
optical and X-ray data sets (see Appendix [A] for details). 



3 SURFACE BRIGHTNESS PROFILES 

In this section we discuss the surface brightness profiles of 
BCGs in the 48 CCCP clusters. We measure the surface 
brightness in circular annuli. The profiles are presented in 
appendix [C] 

We fit circularly symmetric, PSF-convolved r 1 / 4 mod- 
els (deVaucouleurs 1948) to the surface brightness profiles of 
the CCCP BCGs, which yield estimates for the effective ra- 
dius (r e ) and effective surface brightness (fl e ). Details of the 
fitting procedure are given in appendix IB1I Individual fits 
to each system are shown in figures [CTl and [C2l The model 
fits are of varying quality among the CCCP BCGs. Most 
BCGs are well fit, but exhibit deviations from a pure deVau- 
couleurs profile at either small or large radii. We note that 
more complicated models (e.g. Sersic double-r ' ) 

may fit the profiles better but for our purposes the improve- 
ment is negligible. For further discussion on the validity of 
various choices of surface brightness models see Gonzalez et 
al. (2005) and references therein. 

Some of the observed profiles exhibit steep changes in 
the surface brightness over a small range in radius (eg. 
Abell 223a, Abell 521, Abell 697, CL1938+54). These steep 
changes in brightness are caused by imposing a circularly 
symmetric profile on an object which is intrinsically ellipti- 
cal. These sharp variations occur wherever large portions of 
the BCG major or minor axis are masked. For instance, if 
the major axis is heavily masked in a given radial bin, the 
mean surface brightness in that bin is driven towards a value 
that is more representative of surface brightness along the 
minor axis. The impact of this effect is most noticeable for 
those objects which are highly elliptical (CL1938+54) and 
also in the case where opposing sides along the major/minor 
axis are heavily masked in the same radial bin. 

3.1 Kormendy Relation 

The Kormendy Relation (KR) (Kormendy 1977) is an em- 
pirical scaling relation between /i e and r e for elliptical galax- 
ies and has been shown to be a projection of the more gen- 
eral Fundamental Plane (PF) (Djorgovski and Davis 1987, 
Dressier et al. 1987), which describes a 2-dimensional man- 
ifold in (co,{^)e,r e ) space. The KR was first observed to 
hold for local ellipticals and has since been applied to a 
large range of spheroid galaxy types. 

The KR for the CCCP BCGs is shown in figure [1] Sur- 
face brightness fi e is corrected for the effects of cosmological 
dimming and galactic extinction. Transformations to John- 
son V magnitudes are computed by applying an elliptical 
galaxy template spectrum (Coleman et al. 1980). Only those 

2 http : //www. cf ht .hawaii . edu/Instruments/Elixir/ 
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Figure 1. Kormendy Relation (KR) for the CCCP BCGs. The 
data has been corrected for cosmological dimming and converted 
to rest-frame V magnitude. Black squares denote the best fit r 1 / 4 
model parameters. Error bars show the 1 a fitting uncertainties 
but are typically smaller than the symbol size. The solid black line 
shows the mean relation for the CCCP BCGs at z = after cor- 
recting for passive surface brightness evolution. The dashed, black 
line shows the extrapolated KR for local elllipticals (Kormendy 
1977). The circled points show the BCGs that host blue-cores (see 
section [4j . The upper right panel shows the residual in effective 
surface brightness as a function of redshift. 

parameters from the fits to the r' and R data are used as 
their V magnitude conversions are less affected by the as- 
sumption of a particular spectral template compared to the 
g' and B data. Only systems which are unambiguously se- 
lected as the BCG are shown (i.e. selection code=l, see sec- 
tion[S}. The CCCP BCGs are plotted in black squares. The 
residual in /i e from the mean KR (Afj, e ,v) versus redshift is 
shown in the inset of figure [T] If we assume that the offset is 
due to passive evolution we find a mean surface brightness 
evolution of: 

A/^v = ( - l.U ± 0.40)z + (0.28 ± 0.20) (1) 

The relation [TJ is shown in the inset as a solid line and de- 
scribes a slow dimming of the dominant stellar population 
as the BCG ages. This type of passive evolutionary behav- 
ior is expected for massive galaxies that form the majority 
of their stellar content at relatively early cosmic times. If 
we further make the assumption that BCGs are homologous 
this surface brightness evolution equates to a change in the 
mass-to- light ratio M/L with time: 

AlogM/Lv = -0.44 ± 0.16Az (2) 

This result is in agreement with work by van Dokkum et 
al. (1998) who use FP observations of cluster ellipticals to 
obtain a M/L evolution of AlogM/L B = -0.55 ± 0.05Az 
for a ACDM cosmology. 

We apply the evolution correction (equation [1]) to the 
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data and compute the z = KR for BCGs (thin solid black 
line). The scatter in the evolution-corrected data obcg is 
0.28, which is consistent with the scatter in the KR of local 
ellipticals (Kormendy (1977)). The best fit slope of the KR 
for CCCP BCGs is a BC G = 3.44 ± 0.13, which is steeper 
than the values of clbcg = 3.125 ±0.14 reported by Oegerle 
& Hoessel (1991) and a E = 3.02 ± 0.14 found for local el- 
lipticals. This is qualitatively consistent with the result of 
von der Linden et al. (2007) who argue that the difference 
in the KR slope between BCGs and regular ellipticals may 
be explained by a disparity in M/L or a break from homol- 
ogy between these two types of galaxies. Such a difference in 
M/L indicates an increased dark matter fraction in BCGs, 
while a break from homology implies a difference in the way 
these two types of systems are assembled. 

The points which are circled in figure Q] correspond to 
the BCGs that host blue cores (see section [4}. Most of these 
systems have r e < 35 h^g kpc and in the CCCP sample 
they dominate the small r e end of the KR. The majority of 
the blue-core systems lie on the bright side of the best fit 
to the data, but it is important to note that they do not 
lie significantly off the relation. This implies that despite 
the blue emission in these systems, the surface brightness 
fit is more sensitive to the spatially extended, older stellar 
population. We find that removing the blue-core BCGs from 
the KR leads to a further steepening of the slope to ascG = 
3.66 ± 0.15. The slope of the KR for the CCCP BCGs is 
thus in disagreement with that of local ellipticals, regardless 
of whether or not the blue-core systems are included in the 
analysis. 



4 COLOUR PROFILES 

Using our multi-filter photometry we construct colour pro- 
files of the BCGs in our sample, colour profiles trace the 
distribution of stellar populations in these systems, which 
is important for constraining formation and evolution mod- 
els. As shown in figures IC3I and IC4I in the appendix, the 
general trend seen in BCG colours is a gradual decline to- 
wards the blue with increasing radius. This is expected for 
an age/metallicity gradient in massive ellipticals (Vader et 
al. 1988). 

There are several exceptions, however, as 25% (13/53) 
of BCGs in the CCCP show a reversal from this basic behav- 
ior with the colour becoming increasingly blue in the core. 
These galaxies are displaced from their host-cluster red se- 
quence by ~ 0.5 to 1.0 mag in (g'-r 1 ), which is important 
for optical cluster studies that may reject the BCG based 
on colour. The blue cores in these systems cannot be fully 
explained by AGN point-sources because they are at least 
twice the size of the seeing disk. Two of the more extreme ex- 
amples of blue-core BCGs are those in Abell 1835 and Abell 
2204. These systems are hosts of CO emission in their central 
regions (Edge 2001) and are found to host significant recent 
star formation (Crawford et al 1999, Wilman et al. 2006). 
Four other BCGs in the MegaCam sample that show simi- 
lar trends in their colour profiles are: the BCG in RXJ1347 
(ie., the brightest known cluster in X-ray (Schindler et al. 
1995)), the BCG in 3C295, which shows visible signs of 
strong AGN activity in X-ray and radio as well as signs 
of massive amounts of star formation (Thimm et al. 1994, 
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Figure 2. BCG M 1 " 1 vs. BCG/X-ray offset. Blue-core BCGs 
are identified with large blue circles. The BCGs with host clus- 
ter T x < 5 keV are shown with open red triangles. There is an 
obvious tendency for brightest BCGs to lie closest to their host 
cluster's X-ray peak while the faintest BCGs are the furthest. The 
starred symbol represents the BCG in Abell 2261 and is believed 
to host dust enshrouded star formation. The anomalously faint 
blue-core BCGs belong to Abell 2390 and Abell 115N (see text 
for discusssion). 

Harris et al. 2000), the BCG in CL0910 which is identified 
as a hyper- luminous infrared galaxy (HLIRG; Li r > 10 13 L© 
Goto 2005), and the BCG in A115N which is a confirmed 
host of recent star formation (Crawford et al. 1999). Of the 
BCGs in the CFH12K sample Abell 68, Abell 370a, Abell 
383, Abell 2390, MS1455 and MS1512 are classified here as 
blue-core BCGs. 

Table [3] lists star formation rates derived from a litera- 
ture search for a subset of blue-core BCGs. The SFRs range 
from ~ 1 to ~ 100 Mq yr _1 . Although the table indicates 
that the typical SFR of a blue-core BCG is significant, it is 
important to note that it is an order of magnitude less than 
the SFR predicted by the classical cooling flow model. 

4.1 Location within the Cluster 

If the central star formation in the blue-core BCGs arises 
from gas that is cooling out of the ICM, then we expect 
these systems to be located exclusively at the centres of 
relaxed clusters, where the cool material is deposited. To 
investigate the influence of cluster-centric distance on the 
luminosity of BCGs and the presence of blue-cores, we ex- 
amine relationship between BCG total magnitude and the 
projected distance of the BCG from the cluster X-ray peak. 

The X-ray peak positions are obtained using archival 
images from Chandra and XMM-Newton. The total B-band 
magnitude of the BCG {M B ot ) is calculated by integrating 
the measured r' or R surface brightness profile out to r e and 
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Table 3. Star formation rates derived from a literature search 
for a subset of blue-core BCGs 
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Ref. 










3C295 


5.0 


[Oil] & [OIII] 


(1) 


Abell 115N 


0.2 


SSP (spectra) 


(2) 




6.2 


IR 


(3) 


Abell 1835 


125 


IK 


(4) 




40 


Ha 


(4) 




77 


SSP (spectra) 


(2) 




140±40 


SSP (U-R) 


(6) 




123±5 


UV excess 


(7) 


Abell 2204 


14.7 


IR 


(3) 




1.3 


SSP (spectra) 


(2) 


Abell 2390 


5 


IR 


('1) 




5 


Ha 


(4) 




5.4 


SSP (spectra) 


(2) 


CL0910+41 


41±12 


SDSS [0II]ew 


(5) 


MS1455+22 


8.55 


SSP (spectra) 


(2) 




36±4 


UV excess 


(7) 



(1) Thimm et al. 1994, A&A 285, 785 

(2) Crawford et al. 1999 MNRAS, 306, 857 

(3) Quillen et al. 2007, larXiv:0711.TTl8l l 

(4) Egami et al. 2006, ApJ 647, 922 

(5) SDSS DR6 

(6) McNamara et al. 2006, ApJ, 648, 164 

(7) Hicks & Mushotzky 2005, ApJ, 635L, 9 

Star formation rates (SFR) for several blue-core BCGs 
in the CCCP sample. Technique refers to the method used 
to calculate the SFR. For references (3) & (5) the SFR was 
calculated using the indicated measure and the appropriate 
Kennicutt SFR conversion. The rates shown here are signifi- 
cant, ranging from ~ 1 to ~ 100 Mq yr — 1 . It is important to 
note, however, that the SFRs assembled here are an order of 
magnitude smaller than the estimates of the classical cooling 
flow model, which expected SFR ranging from ~ 100 to 
~ 1000 M Q yr" 1 

doubling the resulting flux to account for the light outside 
r e . This total red flux is then converted to rest-frame B 
magnitude, assuming an elliptical galaxy template spectrum 
of Coleman et al. (1980). 

Figure [2] shows the BCG total B-band magnitude 
(Mg ot ) as a function of the projected offset from the host 
cluster X-ray peak (R a)- The BCGs with blue cores are 
identified with large, blue circles. It is immediately clear 
from this figure that all blue-core BCGs have small R a s; 
90% (9/10) of the blue-core BCGs have R oS < 10 kpc. 
This is highly suggestive that the recent star formation in 
these systems is related to core processes within the cluster. 
The system marked by a red star in this figure represents 
the BCG in Abell 2261. From analysis of sub-mm images 
taken with SCUBA on JCMT, Chapman et al. (2002) find 
a significant amount of dust (9.2 x 10 7 h^o Mq) associated 
with the BCG in A2261. Our high-resolution imaging of this 
object reveals a double nucleus feature, possibly created by 
the presence of a dust lane. These observations suggest that 
although this system lacks a blue core, it is either a host of 
dust-enshrouded star formation or that the conditions are 
ideal for the ignition of a starburst in this galaxy (ie. a large 
reservoir of cold gas). 



The three systems marked with large, open, red trian- 
gles are the clusters with T x < 5keV. Because cluster mass 
and X-ray temperature are correlated, we expect these sys- 
tems to be the least massive clusters in our sample. Further- 
more, since BCG luminosity scales with cluster mass (Lin & 
Mohr 2004) these three BCGs are expected to be under- 
luminous compared to their counterparts in more massive 
clusters. Despite the dynamic mass range of the CCCP sam- 
ple, figure [2] shows a general trend of increasing A/g ot with 
decreasing R a g. There is a notable lack of bright BCGs with 
large Rod as well as faint BCGs with small i? ff . This may 
be understood considering the following: the dependence of 
dynamical friction on satellite mass causes the more mas- 
sive BCGs to sink to the centres of their clusters on shorter 
time scales than the less massive BCGs. Furthermore, BCGs 
that are located in cluster centres grow more rapidly than 
those in the outer regions of cluster because they have an 
enhanced probability of trapping and cannibalizing smaller 
galaxies and in cool core clusters they are the beneficiaries 
of the cooling flows. The combination of these effects will 
give rise to the trend in figure [2] 

The obvious outliers from the trend in figure [2] are the 
blue-core BCGs in Abell 2390 at « -21.7 and Abell 

115N at Mg* w -21.5. These systems have best-fit r e val- 
ues of 18 and 15 kpc respectively; quite small compared to 
the rest of the BCG sample. Neither of these galaxies have 
a well developed stellar envelope which could be the rea- 
son they are well fit with such compact surface brightness 
models. The optical images of Abell 2390 show a wealth of 
substructure surrounding the BCG. Bardeau et al. (2007) 
find that the weak lensing signal in this cluster requires an 
elongated mass mass distribution that matches the elonga- 
tion seen in high-resolution X-ray imaging. These results 
support the notion that the BCG in this cluster is young 
and has yet to cannibalize enough substructure to form a 
large cD envelope. Abell 115N is the northern BCG in a 
binary cluster that appears to be in the early stages of a 
merger. Metevier et al. (2000) find that this cluster has an 
anomalously high blue fraction for its redshift. It may be 
that the reason the BCG in Abell 115N lacks a cD envelope 
is related to the age of the remainder of the galaxy popula- 
tion, which may not yet have had time to dynamically relax 
and become cannibalized. 

We perform a Kendall's r test on these data, which re- 
veals that Mg and R s are correlated with a r = 0.26 at 
the 98% confidence level (~ 2.5cr significance). The inter- 
pretation of the trend seen in figure [2] is complicated by the 
possible presence of selection effects, especially at the faint 
BCG/large offset end of the diagram. Some of the clusters 
in the CCCP sample are in the early stages of a merger 
between two or more sub-clusters. These systems naturally 
have large values of R s and do not have well developed cD 
galaxies. Future work involving mass modeling of the clus- 
ters from the combined X-ray, optical, S-Z and weak lensing 
measurements will allow a more detailed investigation into 
these effects. It is clear from examination of the X-ray pro- 
files of the largest R g systems however, that not all of these 
objects are unrelaxed merging sub-clusters and thus the re- 
lation is not completely due to this type of selection bias 
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Figure 3. The L x -T x relation and BCG inner colour gradient for CCCP clusters. A simple measure of the inner colour gradient is 
obtained by performing a linear fit to d{g' — r')/dlog(r) (MegaCam) or d(B — R)/dlog(r) (CFH12K) between at r = 20 ft.^, 1 kpc and the 
FWHM of the limiting PSF. Colour gradients in (B — R) are converted to (g' — r') using the transformations of Fukugita et al. (1996). 
Left: L x -T x relation for the CCCP cluster. Blue circles indicate cluster that contain BCGs with inverted central colour gradient (blue 
cores). Right: colour gradient vs. offset in L x from the mean L x -T x relation at a given T x : Aj /i = Alog(L x /E(z) ergs -1 ) (X-ray 
excess). Note the blue core BCGs are clustered at positive values of X-ray excess. 



5 RELATION TO X-RAY PROPERTIES OF 
THE HOST CLUSTERS 

In this section we compare the optical properties of the 
BCGs in our sample with the X-ray properties of their host 
clusters. The systems in the CCCP sample have been ob- 
served with ASCA and their X-ray temperatures and lumi- 
nosities are taken from Horner (2001). 

There is increasing evidence that the intrinsic scatter 
in the L x -T x relation is physical in origin and that the pro- 
cesses that give rise to the scatter (ie. radiative cooling, star 
formation/stellar winds, AGN) are also responsible for the 
deviation in M-L x as well as SZ correlations (Babul et al. 
2002, McCarthy et al. 2003, McCarthy et al. 2004, Balogh 
et al. 2006). 

McCarthy et al. (2004) demonstrated that models of 
cluster evolution that simultaneously incorporate radiative 
cooling along with a source of preheating reproduce a sig- 
nificant amount of the observed scatter in the cluster L x -T x 
relation. In these models, preheating drives a cluster of a 
given mass towards lower X-ray luminosities. This creates 
an entropy floor in the cluster core which limits the central 
density and therefore L x . Radiative cooling can erase the 
entropy floor and produce dense, bright L x cluster cores. 
These cool core clusters are expected to lie preferentially on 
the high luminosity side of the L x -T x relation. 

An interesting question is whether the same processes 
that are believed to cause the scatter in the L x -T x relation 
also alter the optical properties of BCGs. The CCCP is an 
excellent sample to investigate this with because it spans the 
full observed range of scatter in the cluster L x -T x and M- 



T x relations. To test this scenario we compare the deviation 
of cluster L x from the best fit L x -T x to various properties 
of BCG. The quantity A Lx = Alog(L x /E(z) h^ergs' 1 ), 
with E(z) = (r^Af (1 + ^) 3 + ^a) 1//2 , represents this deviation 
corrected for the change in cosmology and the evolution of 
the mean background density. In what follows this metric is 
used as a first-order proxy for the cooling state of a cluster. 



5.1 Blue-core BCGs and the L x -T x relation 

In section 14.11 we found that the BCGs hosting recent star 
formation are all located near the peak of their host cluster 
X-ray distribution. We now seek to discover the differences, 
if any, between clusters that host blue-core BCGs and those 
that do not. With this aim, we examine the inner colour 
gradients of the BCGs in the CCCP in relation to the cooling 
state proxy A_t x described in section 

The L x -T x relations for the CCCP clusters is shown in 
the left hand panel of figure[3] The data have been corrected 
for any differences in cosmology and scaled by E(z) to com- 
pensate for the evolution in the cosmic background density. 
The mean relation, shown as a dotted line, is determined 
by fitting a power law of the form L x /E(z) = (3T X to the 
clusters with regular BCGs (ie. red cores). The large blue 
circles show those clusters that contain BCGs with positive 
central colour gradients (blue cores); where colour gradient 
is estimated by performing a linear fit to d ^ a ~^ or d ^~f^ 
between r = 20h n l kpc and the FWHM of the limiting PSF. 
Cool-core clusters are known to populate the high L x side of 
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the scatter in the L x -T x diagram, which is where the clusters 
that contain blue-core BCGs are found. 

The right hand panel of figure [3] shows the relationship 
between Al x of the host cluster and the central colour gra- 
dient of the BCG. colour gradients measured in (B — R) are 
converted to (</ — r') using the transformations of Fukugita 
et al. (1996). The small red symbols represent BCGs with 
negative central colour gradients (red cores) and the large 
blue circles represent the BCGs with positive central colour 
gradients (blue cores). As seen in the left hand panel of 
figure [3j all of the BCGs with positive colour gradients (a 
sign of recent star-formation) are found in clusters with the 
highest values of X-ray excess. Moreover, above Ai B > 0.25 
all clusters contain BCGs with positive central colour gra- 
dients. The apparent requirement that clusters have an L x 
excess above the mean L x -T x relation in order to host blue- 
core BCGs implies that the source of the cold gas required 
to fuel these star-forming blue-core systems is linked to the 
processes that give rise to the L x excess. This result there- 
fore points to the cooling ICM as the likely source of cold gas 
replenishment and not accretion from an interactions with 
gas rich substructure. Furthermore, the ubiquity of star- 
forming cores in clusters with sufficiently high Al x suggest 
that heating in these systems (ie. AGN feedback) generally 
tempers the radiative cooling but does not shut it off alto- 
gether. 

5.2 Scatter in L x -T x and BCG/X-ray peak offset 

The results discussed in the previous sections indicate that 
the presence of blue-cores in BCGs is linked to the BCG's 
position as well as the X-ray excess of the cluster. This 
is consistent with the hypothesis that cooling in the ICM 
is providing the cold gas necessary to fuel the recent star 
formation. This star formation is one mechanism by which 
BCGs can acquire stellar mass, but they also grow via merg- 
ers and accretion. Through the cluster density profile these 
latter two processes also depend on the BCGs cluster-centric 
location, ft is therefore interesting to examine to what de- 
gree the sources of scatter in the cluster L x -T x relation affect 
the position of the BCG within the cluster or vice-versa. 

In figure U we show the dependance of cluster A l x on 
i? ff. The symbols in this figure are the same as those in fig- 
ure[2] Strikingly, all of the clusters with values of Ai B > 0.25 
have BCGs that are within ~ Wh? kpc of their cluster X- 
ray peaks. Furthermore, virtually every system above this 
X-ray threshold host a star-forming BCG. This is expected 
if the cold gas that is fueling the star formation in these sys- 
tems originates through cooling of the ICM. Abell 2261 is the 
only cluster with a large value of Al b and a small value of 
R g that does not contain a blue-core BCG. As discussed in 
section |4~T1 this BCG has a large dust mass 9.2 x 10 7 Mq 
(Chapman et al. 2002). Bonamente et al. (2006) show that 
A2261 has a central cooling time t coo i < 0.5t hubb i e , as well as 
a central temperature gradient, indicating that it has an ac- 
tively cooling core. The existence of a system such as Abell 
2261, where large amounts of cold gas remains unconverted 
into stars indicates that there may be additional considera- 
tions required for understanding star formation in actively 
cooling systems. Perhaps there is a time delay between the 
build-up of fresh cold gas and the onset of star formation. 
It may also be possible that an interaction with a merging 
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Figure 4. A^ x vs. BCG/X-ray offset. Clusters with blue-core 
BCGs are identified with large blue circles. Strikingly, all of the 
BCGs in systems with Ar, > 0.25 are located very close to 
the central X-ray peak and are hosts of recent star formation. 
The remainder of the systems are found at a wide range in R g . 
Only a single system (Abell 2261 - starred symbol) in the small 
offset/high A^ x group does not contain a blue-core BCG; we find 
evidence in the literature, however, that suggests this system is 
host to dust enshrouded star formation 



substructure is required in order to trigger the star forma- 
tion; or that star formation is present but hidden by the 
dust. Regardless of these details, figure [4] demonstrates that 
the clusters with large positive residuals from the L x -T x all 
have their BCGs relatively well aligned with the cluster X- 
ray peaks; indicating an advanced state of dynamical relax- 
ation. 



5.3 Effect on BCG Luminosity 

Recent observations (Crawford et al. 1999, Wilman et al. 
2006) show that some cool-core clusters are hosts of mod- 
erately high star formation rates in their centrally located 
BCGs (eg., Abell 1835: SFR ~ lOOMeyr -1 ; Abell 2204: 
SFR ~ IMoyr -1 ; see table [3}. The radiative losses in 
these systems are therefore not fully balanced by the heat- 
ing sources; instead the energy budget is such that some of 
the gas is cooling enough to form stars. Figure [3] illustrates 
that the BCGs which are forming stars are exclusively lo- 
cated in clusters on the high L x side of the scatter in the 
L x -T x relation, which is the region of the diagram that is 
populated by cool-core clusters. In order to further under- 
stand and quantify the importance of the cooling state of 
a cluster in contributing to the stellar content of BCGs, we 
investigate in this section how BCG luminosity varies with 
the X-ray excess of the host. 

Figure [5] shows BCG Mg l versus A Lx . The point size 
is scaled according to projected offset of the BCG centre 
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Figure 5. BCG total B-band magnitude (M^) versus A Lm of 
the host cluster. Symbols increase in size according to the BCG 
projected offset from the X-ray peak (Raff). Blue core BCGs are 
encircled in blue. Large red triangles identify the clusters with 
T x < 5 keV. These systems are expected to host BCGs that are 
less massive than those in the rest of the sample (see text). There 
is a slight trend of brighter BCGs being found in clusters with an 
X-ray excess compared with the mean L x -Tx relation. Many of 
the outliers from the trend have large offsets with respect to the 
X-ray peak. Small open symbols represent the systems for which 
no X-ray data is available and so no offset is determined. 



from the peak of the X-ray distribution. X-ray peak posi- 
tions are obtained from either Chandra or XMM-Newton 
observations. The bin limits in X-ray offset are 10, 25, 50 
and 100 h^Q kpc. The points marking clusters with blue 
cores are surrounded with blue circles. The clusters lacking 
high resolution X-ray data are plotted with open symbols. 
The few clusters from the archival data sample with T x < 5 
keV are marked with large, open, red triangles. 

Examination of figure [5] reveals a trend of increasing 
BCG luminosity with increasing host cluster deviation from 
the mean L x -T x relation. This trend appears to be some- 
what weak but if the systems with T x < 5keV (ie. those 
shown with large, open, red triangles) are ignored the trend 
becomes more apparent. Furthermore, many of the outliers 
are BCGs which lie at large offsets from their host cluster's 
X-ray peak, which is expected if the L x -T x deviations arise 
from processes in the cluster core. Particularly striking is the 
absence of both faint BCGs in clusters with an L x excess and 
bright BCGs in clusters with an L x deficit. A Kendall's r test 
of this correlation including all data points shown in figure 
[5] yields r = —0.25 with a confidence level of 98% (~ 2.5cr 
significance). The trend displays a full magnitude difference 
for BCGs across the range in observed variation from the 
cluster mean L x -T x relation, implying that the source(s) of 
scatter in the cluster L x -T x relation contributes significantly 
to the growth of BCGs. It is interesting that while only the 



clusters with the highest values of Al x display the observ- 
able signatures of active star formation (see figure [3j , the 
trend in BCG luminosity spans the entire range of observed 
X-ray excess. This suggests that another important mecha- 
nism in addition to recent star formation is contributing to 
the stellar mass differentiation of BCGs. As expected, the 
blue core systems inhabit mainly the upper right portion 
of the diagram. The recent star formation in these systems 
may explain their high average luminosity with respect to 
the rest of the sample. Two exceptions are the low luminos- 
ity BCGs in Abell 115N and Abell 2390, which are believed 
to be young BCGs that have yet to cannibalize a large por- 
tion of the galaxy population in their respective clusters (see 
section 14.11 for discussion) . 

We examine the relationship between Mjj ot and host 
cluster L x as well as Mjj ' and host cluster T x . We find no 
significant evidence for a correlation between either of these 
two pairs of variables. 

The existence of the trend seen in figure [S] (M^ vs. 
A_l x ) suggests that whatever is responsible for producing the 
scatter in the L x -T x relation is also assisting in the growth 
of stellar content of BCGs. We have argued that cooling is 
the primary cause. Alternatively, scatter may be introduced 
into the L x -T x relation by mergers with massive substruc- 
tures as well as by intrinsic variations in the shapes of the 
underlying dark matter halo density profiles from one halo 
to the next. Recent papers by Rowley et al. (2004), Kay et 
al. (2006) and Poole et al. (2007) show that although merg- 
ing in clusters is capable of generating some scatter in the 
X-ray scaling relations, this scatter is short lived and does 
not account for the degree of scatter seen in the observa- 
tions. In addition, work by Balogh et al. (2006) shows that 
the variation in dark matter density profiles also introduces 
some scatter but that the resulting scatter is not sufficient 
to reconcile the observations. These studies therefore add 
support to the interpretation that the variation in BCG lu- 
minosity is directly linked to the balance between heating 
and cooling in the ICM. 



6 CONCLUSIONS 

We have conducted a detailed comparison of the optical 
properties of BCGs with their host cluster X-ray proper- 
ties in 48 intermediate redshift systems using a sample that 
spans the full observed range in cluster morphology. Our 
analysis reveals a close connection between the evolution of 
BCGs and that of their host clusters. The results are dis- 
cussed below. 

The Kormedy Relation defined by the CCCP BCGs is 
found to have a steeper slope than that of regular elliptical 
galaxies (ancG = 3.44 ± 0.13 versus cle = 3.02 ± 0.14). 
This difference in slope may be influenced by the sample 
selection and the use of circular models; but it suggests that 
the formation processes of the CCCP BCGs differ from their 
regular elliptical counterparts. Von der Linden et al. (2007) 
propose that the observed difference between the slope of the 
Fundamental Plane of BCGs and that of regular ellipticals 
may be a result of an enhanced dark matter fraction due to 
the typical position of the BCG near or inside the cluster 
core. Long-slit spectral observations of BCGs in the CCCP 
sample will allow us to better determine the cause of these 
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differences and examine how baryonic feedback in clusters 
affects the scaling relations for galaxies. 

The majority of BCGs in the sample have shallow colour 
profiles that become bluer with increasing radius. There are 
several BCGs (25%) that deviate from this simple behavior; 
and show blue cores instead. This colour peculiarity results 
in large deviations from the red sequence (up to ~ 1.0 mag 
in (g'-r')) and affects optical cluster studies that use colour 
to determine cluster membership. The radial extent of the 
blue-cores shows that they cannot be explained by point 
sources alone. This leads us to conclude that they are re- 
gions of active star-formation. Several of the blue-core BCGs 
have measured star formation rates in the literature (see ta- 
ble [3]). Dust and molecular gas has also been detected in 
many of these systems (Edge 2001, Edge et al. 2002, Egami 
et al. 2006). These recent observations reveal that the evo- 
lution of BCGs differs from the standard picture for massive 
non-BCG ellipticals. The presence of extended star forming 
regions in BCGs requires that the cold gas reservoir is being 
replenished in these systems. Cooling flows in clusters with 
blue-core BCGs provide a natural mechanism to fuel this 
recent star formation. 

One of the predictions of cool core cluster model is that 
the cold gas that condenses out from the ICM is deposited 
in the cluster centre. We find that the star-forming systems 
are all located at small projected offset from the peak of the 
cluster X-ray emission. Having small i? ff appears to be a 
requirement for BCGs to host a blue core; which is consistent 
with the cool core model prediction. 

The competition between baryonic feedback effects, 
such as radiative cooling and heating from stellar 
winds/AGN, is believed to be the source of the intrinsic 
scatter in the cluster L x -T x relation. For a given T x , ra- 
diative cooling increases L x , while heating reduces L x . We 
find that all of the blue-core BCGs are located in clusters 
with large values of X-ray excess; suggesting that cooling is 
dominant in these systems. The observation that all systems 
above Al x > 0.25 contain star- forming BCGs implies that 
all cool core clusters are producing stars in their cores. This 
is somewhat surprising and contrary to the current conven- 
tional wisdom. Furthermore, the lack of star-forming BCGs 
at negative values of X-ray excess is difficult to explain with 
models that fail to take into account the cooling state of the 
host cluster. 

Examination of the relationship between X-ray excess 
and the BCG position within the cluster reveals that all 
of the clusters with > 0.25 contain BCGs (with blue 

cores) that are within ~ 10 kpc from the peak in the X-ray 
emission. This suggests that the conditions responsible for 
the creation of cool-core clusters also lead to the formation 
of a centrally located BCG; both of which are required for 
the presence of a star- forming core. 

The absolute B-bwad magnitude of BCGs is only weakly 
dependent upon the X-ray excess of the host cluster, as well 
as the BCG offset from the X-ray peak. This is not surpris- 
ing given the variety of mechanisms by which BCG acquire 
stellar mass (cannibalism, harassment, star formation) and 
the range in formation times expected from models of hier- 
archical merging. 

These observations support a picture of linked BCG and 
cluster evolution. In this scenario the clusters which contain 
the most dynamically relaxed BCGs are the ones that form 



cool cores and host central star formation. The prototypical 
example of this behavior may be the cluster MS1455+22. 
The optical imaging shows that it is dominated by the cen- 
tral cD galaxy (a blue-core BCG that is aligned with the 
X-ray peak) and contains very few bright satellites galaxies. 
MS 1455+22 appears to be in a significantly more dynami- 
cally evolved state than the majority of the clusters in the 
sample. A similar conclusion is reached by Yee and Elling- 
son (2003) to explain its anomalous position in the velocity 
dispersion-richness plane. MS1455+22 has the highest X-ray 
excess in the CCCP sample; indicating that it hosts a promi- 
nent cool core. We also note that the X-ray imaging shows 
that the ICM in this cluster is highly relaxed. 

Our results are consistent with those of Edwards et al. 
(2007) who find that ~ 70% of BCGs that are within 50 
kpc from their host cluster X-ray peak are line emitting. 
Furthermore, our results are consistent with studies that 
find that the recent star formation in BCGs is associated 
with the ICM cooling rate (Johnstone et al. 1987, Heckman 
et al. 1989, Cardiel et al. 1998). However, by covering the full 
observed range of cluster morphologies, the CCCP sample 
allows us to strengthen these earlier claims and show that 
these phenomena are truly associated with the properties of 
the ICM. 

Taken together, these data have several important im- 
plications. Based on X-ray temperature and luminosity 
alone we can predict that if A_l x > 0.25 the BCG will 
be coincident with the cluster X-ray peak and host a blue 
core. Furthermore, that blue cores are a ubiquitous feature 
of BCGs in cool core clusters (with sufficient X-ray excess) 
indicates that AGN and other heating mechanisms gener- 
ally temper but do not shut off cooling altogether. The data 
also show that BCG photometry alone can be used to pref- 
erentially select cool-core clusters that are prime targets for 
cosmological studies that focus on the cluster gas fraction 

The results presented here are interesting and deserve 
further investigation. Future work requires spatially resolved 
long-slit or IFU optical spectra of the CCCP BCGs. Such 
data would enable one to disentangle the stellar populations 
in these systems and constrain the ages of the starbursts 
in the blue core systems. Furthermore, because the work 
presented here is based solely on photometry we cannot rule 
out some contribution from an AGN component to the blue 
core. Spatially resolved optical and near IR spectra would 
enable one to quantify the contribution of AGN vs. starburst 
as a function of position in these galaxies. 
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APPENDIX A: BCG SELECTION 

The task of identifying the BCG in each cluster is critical 
to the success of this study and for the majority of systems 
(ie., those containing cD galaxies) this task is straightfor- 
ward. These conspicuous galaxies have large, extended stel- 
lar envelopes and are immediately identified as the brightest 
cluster member. It is the clusters which do not contain obvi- 
ous cD galaxies that pose a potential problem. Each of the 
BCGs are given a selection code number from 1 to 3 corre- 
sponding to how easily they are identified, with 1 being an 
unambiguous selection and 3 a highly ambiguous selection. 
Selection codes are listed in tables [1] and [2] 

Al Non-trivial BCG Selections 

There are clusters for which two or even three galaxies are 
potential BCG candidates. In order to select one of these 
candidates over another the additional qualification must 
be made that the galaxy of interest is the candidate which 
is most likely to have been influenced by the cluster's grav- 
itational potential throughout its evolutionary history. This 
additional qualification allows us to take advantage of the 
outstanding high-resolution X-ray observations available in 
the Chandra and XMM-Newton archives. Using these data 
we select the candidate that is located closest to the X-ray 
centroid of its host cluster as the BCG. For a relaxed cluster, 
the X-ray centroid is expected to trace the deepest point in 
the cluster potential. In the following discussion we describe 
the systems with multiple BCG candidates and outline the 
reasoning used in the final BCG selection. 

• Abell 370: The X-ray emission from this cluster is peanut 
shaped with two lobes connected by a bridge. Each lobe has 
an associated BCG candidate. This system is likely a pair 
of sub-clusters in the early stages of a merger and thus we 
treat both candidates as individual BCGs with the north 
(N) and south (S) BCGs labeled 'a' and 'b' respectively. 

• Abell 520: This cluster has a highly disturbed X-ray mor- 
phology indicative of a major merger. There are two BCG 
candidates lying on the outskirts of the cluster X-ray con- 
tours. The NE candidate is chosen as the BCG over the SW 
candidate because it lies closer to the centroid of the overall 
X-ray distribution. 

• Abell 959 contains several BCG candidates. There are two 
galaxies in particular that lie at small projected distances 
from the X-ray peak. Of these two we choose the N candi- 
date as the BCG because it lies slightly closer to the X-ray 
centroid. We note that this selection is somewhat uncertain 
because the X-ray image that is used to arrive at this selec- 
tion is from Rosat. 

• Abell 1234: There are two potential BCG candidates in 
this clusters. They are identified as the NW candidate and 
the SE candidate. There are no X-ray data for this cluster 
in either the XMM-Newton or Chandra archives. We select 



the NW candidate as the BCG because it is more centrally 
located within the cluster and it appears to have a more 
extended stellar envelope in comparison with other candi- 
dates. 

• Abell 1246: This cluster contains two potential BCG can- 
didates. These object are overlapping on the sky and distin- 
guished as the N candidate and S candidate. The S candidate 
has a higher peak surface brightness. The N candidate how- 
ever has a position angle that is more aligned with the over- 
all extended profile of these two components. The extended 
profile also seems to be more centred on the N candidate. 
We choose the N candidate as the BCG. 

• Abell 1914: This cluster has two BCG candidates with one 
in the NE and the other in the SW. The SW candidate has 
a brighter peak surface brightness than the NE candidate 
but it lies further from the X-ray centroid by a significant 
margin. We choose the NE candidate as the BCG. 

• Abell 2111: This cluster has a pair of BCG candidates. 
We select the N candidate as the BCG over the S candidate 
because it is closer to the X-ray peak and it has a much 
more extended low surface brightness envelope. 

• CL0910+41: This cluster has a very bright BCG that 
sits near the centre of the optical light distribution 
(SDSSJ091345. 5+405628). This galaxy is unlike many other 
BCGs in that it is compact, extremely blue in colour and 
contains very little evidence of an extended stellar envelope. 
These differences raise the concern that this galaxy may 
be a foreground source that projects onto the cluster. Goto 
(2005) identifies this galaxy as a Hyper Luminous InfraRed 
Galaxy (HLIRG) with an infrared luminosity Lt r > 10 13 L© 
at z — 0.442 which makes this galaxy a confirmed cluster 
member and validates our BCG selection in this case. The 
extraordinary properties of this galaxy as a BCG make it 
particularly interesting. 

• MS0906-11: The galaxy chosen as the BCG in this cluster 
was selected based on its small offset from the central X-ray 
peak as well as the flat shape of its surface brightness profile 
compared to other potential BCG candidates in this cluster. 

• MS1621+26: The X-ray image shows a disturbed morphol- 
ogy containing an elongated central peak with a plume of 
emission extending from the central peak to the north. Of 
the two BCG candidates (N and S) we select the southern- 
most candidate because it lies closer in projection to the 
central peak while the northernmost candidate appears to 
be more associated with the plume of emission. It is also in- 
teresting to note that the selected BCG has a double peaked 
core visible in the R image. This system appears to be in an 
intermediate stage of a merger between two sub-clusters. 



A2 Clusters with strong BCG Ambiguities 

After selecting BCGs in clusters for which the selection is 
not immediately obvious there still remains a group of clus- 
ters in which multiple BCG candidates are similar enough 
according to our selection criteria that it is not possible to 
reliably label a single galaxy as the dominant BCG. In the 
following we briefly describe the systems for which this is 
the case. 
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• Abell 851 contains several candidate BCGs. We narrow the 
selection to a closely spaced group of three bright galaxies 
that lie nearest the centre of X-ray emission in an XMM- 
Newton 24 ks exposure. Of these three remaining candidates 
we select the southernmost galaxy as the BCG because it 
appears to have the largest extended stellar envelope. 

• CL0024+16: This cluster has a complex core with four 
bright, red galaxies spaced with their surface brightness pro- 
files significantly overlapping such that they appear to share 
a common diffuse halo. It is likely that these galaxies are 
in the processes of merging to form a cD galaxy. A 40ks 
Chandra image of this field shows a mostly relaxed X-ray 
surface brightness profile which is elongated in the direction 
of the four bright, red, central galaxies. The easternmost 
of the four red galaxies is coincident with the X-ray peak 
however the active dynamical state of the system makes the 
selection of the BCG based on the X-ray distribution less 
reliable because of possible separation between the gas and 
the red galaxies. The combined difficulty in BCG selection 
and modeling of the overlapping red galaxy surface bright- 
ness profiles cause this cluster to be flagged and removed 
from some of the analyses. 

• MS0016+16: This cluster contains three potential 
BCG candidates which in projection are seen to lie along 
a line from the northeast to southwest. The two most north- 
east candidates lie closest to the central X-ray peak of the 
cluster but are both at the same projected clustercentric 
distance. It is likely that these three systems will eventually 
merge to form the BCG but in the current dynamical state 
and viewing angle its is not possible to reliably select a BCG 
from this cluster. 

• MS1231+15: There are two obvious BCG candidates 
in this cluster which appear to be in the intermediate stages 
of merging together. They share a common halo of diffuse 
emission that is elongated along the axis connecting the two 
galaxies. There is no X-ray image of this cluster in either 
the Chandra or XMM archives and thus it is impossible to 
select a BCG based on the shape of the X-ray distribution of 
this cluster. These two candidates are labeled MS1231a and 
MS1231b (north and south respectively) but are left out of 
some of the analyses because of the large ambiguity in BCG 
selection and the lack of X-ray imaging. 



APPENDIX B: ANALYSIS 

We remove the image background using a standard 2-d poly- 
nomial model. While a more complex background subtrac- 
tion method would improve the depth of our surface bright- 
ness measurement, such an improvement does not signifi- 
cantly affect the values of the BCG fit parameters. 

For each cluster image we generate a mask which is used 
to reject pixels from inclusion in the surface brightness mea- 
surement of the BCG. The mask is based on a combination 
of 3 layers and is designed to mask everything but the pixels 
that are unambiguously identified with BCG or background 
light. 

The first layer of the mask is derived from a SExtrac- 
tor segmentation image. The initial choice of extraction pa- 
rameters are based on those used in Simard et al. (2002). 



The segmentation image produced is then modified by ex- 
tending each masked region into the background by 3 pixels 
on all sides in order to reduce the contamination from the 
unmasked portion of the detected sources lying below the 
detection threshold. 

The second layer of masking is generated by hand and 
is used to exclude regions of the image that are obviously 
contaminated by sources and defects that are not detected 
by automated source extraction routines such as SExtrac- 
tor (Bertin & Arnouts 1996). Most of the image stack is 
constructed from at least 4 separate exposures however, due 
to small gaps between adjacent chips on the mosaic detec- 
tor there are sections of the stack which contain only 3 or 
fewer exposures depending on the dither pattern used. With 
fewer than 4 exposures it is difficult to determine a robust 
median value when stacking and as a consequence there are 
more cosmic rays present as well as abrupt changes in the 
signal-to-noise properties across these regions. To avoid com- 
plications in surface brightness fitting these regions are also 
masked in this layer. 

Finally, the central regions around the BCG are masked 
by hand out to a radius of 45" for each target BCG. Each im- 
age is carefully examined and all visible sources are masked 
out to the radius at which they become indistinguishable 
from the background level. The contrast is adjusted slightly 
and the process is iterated to fainter and fainter surface 
brightness. The resulting mask portion is used as the third 
mask layer and replaces the isophotal layer (2 nd layer) within 
45" of the BCG centroid. This third layer of masking is re- 
quired to ensure that galaxies in the highly crowded cluster 
cores do not contaminate the fit. Isophotal masking can- 
not account for the shapes of the core galaxies and so it is 
more effective to mask them by hand as the eye can trace the 
falling surface brightness profiles out to large distances quite 
accurately. We have also experimented with modeling and 
subtracting the core galaxies prior to BCG surface bright- 
ness measurement. As Patel et al. (2006), however, we find 
that this technique does not improve the surface brightness 
measurements over those obtained by simple masking. 

Bl BCG surface brightness modeling 

For each BCG surface brightness profile a PSF-convolved, 
circularly-symmetric, single-component r 1 ^ 4 deVaucouleurs 
model is fit with the following functional form: 

I(r) = J e cxp (-7.67((r/r e ) 1/4 - 1)) + I bgd (Bl) 

with I e and r e as the effective intensity and effective radius 
respectively. The term I bg d is a constant background inten- 
sity term which is also included as a fit parameter. To deter- 
mine the PSF we measure the surface brightness profile of 
an unsaturated star close to the BCG. The total integrated 
distribution is then normalized to unity in order to conserve 
total flux when applying the PSF as a convolution response 
function in BCG modeling. We fit a moffat profile (Moffat 
1969) to the PSF star with the following form: 

I(r) = I {l + (r/a) 2 ]-P + I bgd (B2) 

where I bg d is a term added to take into account the local 
background level and (5 is left as a free parameter. 

Estimates for the error in surface brightness for each 
radial bin are obtained by bootstrapping the pixel values in 
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Figure Bl. Marginalized \ 2 contours for surface brightness fit parameters on r' (left) and g' (right) images of Abell 2259. The contours 
are shaded according to confidence level: lcr in red, 2cr in orange and 3<r in yellow. The most notable degeneracy is between r e and I e . 
We note that the dependence of the fit parameters on wavelength is due to the presence of a colour gradient in this system. 



that bin. BCG source pixels are randomly selected in each 
radial bin, allowing a given pixel to be selected multiple 
times. This has the advantage that we do not need to make 
any apriori assumption of the functional form of the noise 
(ie. poisson noise), or to know the precise values of exposure 
times, gain or read-noise. Bootstrapping not only accounts 
for the random noise in the data, it also incorporates some 
of the contribution from non-uniform, systematic variations 
(bright stars, variability in galactic dust absorption, crowd- 
ing effects of nearby galaxy/group light, etc.) across the clus- 
ter image. We note, however, that bootstrapping estimates 
misrepresent the true error in radial bins that have only a 
few pixels available to randomly sample. 

The BCG fitting procedure is restricted to data points 
at radii larger than twice the FWHM of the derived moffat 
profile PSF (r/it > 2r/ u ,; lm ). Ignoring the central bins helps 
to ensure that any small error in the PSF determination will 
not seriously affect the BCG fit. An upper radial cutoff to the 
fitting region is also applied. The fit is initially restricted to 
those data points that lie at radii less than 500/ifQ 1 kpc from 
the BCG centre. This radius is chosen somewhat arbitrarily 
but it reflects an effort to isolate the BCG component within 
the overall cluster diffuse light profile. The PSF-convolved 
deVaucouleurs + background model is then fit to the data 
lying within these two limits as described above. A second 
iteration of the fit is then performed with the maximum fit 
radius set to l.5r e j n i tia i and an extra term is added to \ 2 ■ 
The second fit iteration uses a \ 2 calculated as follows: 

2 2 / Ibgd.new ibgd,initial \ /"Do\ 

Xnew Xinitial I I I ( -Do) 

\ &bgd,initial J 

where the subscripts initial and new refer to the first and 
second fit iteration parameters respectively. The extra back- 
ground term in Xnew acts as a penalty to any major change 



in the background parameter and ensures that the back- 
ground is primarily determined/constrained by the larger 
radius fit range. The reduction in the fit range of the second 
fit iteration helps to further isolate the BCG signal with an 
upper radial cutoff that is motivated by the shape of the 
galaxy itself. 

Examples of the marginalized \ 2 contours are shown 
in figure [5T] for the fits to the Abell 2259 BCG surface 
brightness profile. It is immediately seen that I e and r e are 
the most degenerate parameter pair, highlighting a common 
trait of all of BCG surface brightness fits in the CCCP sam- 
ple. By leaving the background as a free parameter in the 
fitting routine, even though the image has been background 
subtracted, this technique shows how the remaining struc- 
tural properties of the BCG are affected by a possible error 
in background. By inspection of the \ 2 contours shown in 
figures [Bll it can be seen that a 3<r error in the background 
can affect the derived r' effective radius r e at the 3-4 % level 
for Abell 2259. The marginalized \ 2 surface diagrams are 
also useful for analyzing the stability of the fitting routine 
for a given BCG. In summary, these data indicate that our 
BCG fitting procedure determines the structural parameters 
with high precession (within a few %) and that the effects 
of error in the background subtraction are negligible. 



APPENDIX C: OPTICAL PHOTOMETRY 

Figure [Cll shows the surface brightness profiles in the g' and 
r' bands for all BCGs that were observed with MegaCam 
while figure IC2I shows the surface brightness profiles of all 
BCGs that were observed with CFH12K. The measured data 
are shown by the points, solid lines correspond to best fit 
PSF-convolved r 1 ^ 4 models as described in the above section 
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Figure CI. Surface brightnesses in g' (blue) and r' (red) for the MegaCam BCGs in the CCCP sample. Points correspond to measured 
data while solid and dotted lines correspond to fitted r 1 / 4 models with and without PSF-convolution respectively. Vertical dashed lines 
denote the rf w hm seeing limit of the derived moffat PSF. Residuals are shown along the bottom with dotted lines showing the difference 
between fits with and without PSF-convolution and solid lines with error bars showing the difference between PSF-convolved models 
and the data. The top axis shows the corresponding radii in arcseconds 
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Figure C2. Surface brightnesses in B (blue) and R (red) for CFH12K BCGs in the CCCP sample. Points correspond to measured 
data while solid and dotted lines correspond to fitted r 1 / 4 models with and without PSF-convolution respectively. Vertical dashed lines 
denote the rf w hm seeing limit of the derived moffat PSF. Residuals are shown along the bottom with dotted lines showing the difference 
between fits with and without PSF-convolution and solid lines with error bars showing the difference between PSF-convolved models 
and the data. The top axis shows the corresponding radii in arcseconds 
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Figure C3. colour profiles in (g'-r') for BCGs observed with MegaCam in the CCCP sample. Points with error bars indicate measure- 
ments. Solid lines show the difference between best fit g' and r' r 1 / 4 models. Vertical dashed lines show the seeing limits rt w h,m f° r the 
g' data (blue) and r' data (red). The top axis shows the corresponding radii in arcseconds 
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Figure C4. colour profiles in (B-R) for BCGs observed with CFH12K in the CCCP sample. Points with error bars indicate measurements. 
Solid lines show the difference between best fit B and R r 1 / 4 models. Vertical dashed lines show the seeing limits rf w f lm for the B data 
(blue) and R data (red). The top axis shows the corresponding radii in arcseconds 
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and dotted lines show similar fits without PSF convolution. 
Surface brightness data is colour coded by wavelength band 
with r and R data shown in red and g and B data shown 
in blue. In all cases the BCGs are brighter in the red filter 
than in the blue filter with a typical offset on the order of 
one magnitude at all radii. Vertical dashed lines show the 
moffat PSF fwhm radius {rj w hm)- Residuals are shown along 
the bottom of each plot. In each sub-panel the dotted line 
corresponds to the difference between the PSF-convolved 
and non-PSF-convolved fits, and the solid lines with error 
bars show the difference between the PSF-convolved mod- 
els and the measured data. Note that the B-band data for 
MS00I6 has been omitted because it was taken under non- 
photometric conditions. See section [3] for discussion. 

Figure IC3I shows the colour profiles in {g'-r') for the 
MegaCam BCGs, while figure [C4] shows the colour profiles in 
(B-R) for the CFHI2K BCGs. As before, the vertical dashed 
lines denote the rf m hm of the PSFs in r' or R (red) and g' or 
B (blue) . Points with error bars indicate the measured data 
and solid lines show the difference between the red and blue 
r 1 / 4 model fits to the surface brightness distribution. Note 
that even though this is not a direct fit to the colour profile 
the difference between the individual fits to the data frames 
match the colour data well in most cases. See section 2] for 
discussion. 



APPENDIX D: QUANTIFYING SYSTEMATICS 

The measurement of surface brightness profiles of BCGs re- 
quires a detailed understanding of the various sources of 
systematic effects that may contaminate the signal. Previ- 
ous BCG and ICL studies have employed a wide variety of 
tests that aim to quantify these sources of contamination 
(Gonzalez et al. 2005, Krick et al. 2005, Patel et al. 2006). 
The following sections describe the tests that are used in this 
paper to investigate the contributions of systematics effects 
and the steps taken to minimize them. 

Dl Mask Effects 

The measurement of BCG surface brightness is complicated 
by the effects of crowding. The shallow surface brightness 
profiles of other luminous galaxies in the field often overlap 
with that of the BCG. The method used to create the mask 
therefore systematically affects the measured profile. The 
most common mask generation method used in the litera- 
ture rejects pixels around non-BCG galaxies that lie inside 
a particular isophote fii SO (Lin and Mohr 2004, Feldmeier 
et al. 2004, Zibetti et al. 2005, Krick et al. 2005). The ad- 
vantage of this technique is that the methodology is clearly 
defined and reproducible, the disadvantage however, is that 
it does not take into account the profile shapes. The mask 
generation method used here is not limited to a particular 
isophotal cutoff and faint shallow wings of BCG neighbors 
are masked extensively. In this section we discuss the tests 
used to evaluate the influence of mask effects on surface 
brightness. 

The test described here is designed to look for signal 
contamination by comparing the profiles in different angu- 
lar sections of the BCG. The image is divided along the 
BCG major and minor axes into four quadrants. Under the 



assumption of zero contamination, the profiles measured in 
each quadrant should be identical. The differences between 
the profiles in different quadrants is a measure of the asym- 
metric component of the contamination signal. There are 
other causes of variation in the quadrant profiles such as 
asymmetric masking and position angle twists. For instance, 
if an individual quadrant is heavily masked near the minor 
axis then the measured profile in that quadrant is more rep- 
resentative of the surface brightness along the major axis 
and is brighter than the ideal unmasked, zero contamina- 
tion case. In this way, slight variation in the masking from 
one quadrant to another introduces scatter among the four 
quadrant profiles. This test is applied to all BCGs and the 
presence of large scatter in the quadrant profiles is used to 
determine the radial limit to which the BCG surface bright- 
ness data are reliable in each of the images (rdat,max)- 

The plots shown in figure ID II show the result of the 
masking incompleteness tests for the BCGs in the clusters 
Abell 2261 and Abell 2537. The profiles in the four quadrants 
are shown as different coloured curves. Dotted lines indicate 
the radial limit (rdat,max) to which the the BCG profile is 
minimally affected by asymmetric mask incompleteness. The 
profile deviations seen in the r' image of Abell 2537 just 
inside rdat.max are due to a position angle twist of the BCG. 
This effect is symmetric in two of the four quadrants and 
anti-symmetric in the other two. 



D2 Large Scale Background and Flatfield Error 

Another major source of systematic error in BCG/ICL stud- 
ies comes about from large scale flatfielding error and large 
scale variations in the background. Feldmeier et al. (2002) 
and Patel et al. (2006) present methods that can be used in 
order to quantify some of the properties of background vari- 
ation on a variety of scales. Similar methods are used here 
in order to quantify the image background variability. Each 
image is portioned into a rectangular grid. The background 
statistics are determined within each grid square using the 
same mask that is used to measure BCG surface brightness. 
This process is repeated with variety of grid sizes in order 
to examine background variation on different scales. The re- 
sults of this procedure reveal the presence of structure in the 
background as a function of distance from the BCG centre 
and size of the grid. 

To quantify the effects of large-scale structure we exam- 
ine the statistics of 10 independent surface brightness pro- 
files that are centreed on various parts of the image back- 
ground. From the analysis of several cluster images, we con- 
clude that the variance in the background signal begins to 
dominate the surface brightness measurements at a typi- 
cal limit of fi = 28.5 mag/(") 2 but that brighter than this 
limit, large-scale error does not significantly contribute to 
the BCG signal. 

We test for the possibility that some of the BCG/ICL 
might be removed in the background subtraction process 
due to its large extent. Examining a range of mesh sizes 
used for background determination, we conclude that this is 
not an important effect for the results presented here (5% 
difference in r e in the most extreme case). 
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Figure Dl. Testing for systematic error introduced by mask effects in Abell 2261 and A2537. The BCGs in these clusters are divided along 
their major and minor axes into four quadrants. Each of the different coloured curves in these diagrams represent the surface brightness 
measured in a single quadrant. These diagrams are used to determine the radial range in which the flux can be reliably assigned to the 
BCG (see discussion). Large differences in the profiles from the four quadrants indicate asymmetric sources of contaminating flux. The 
panels (a) and (b) show respectively the r' and g' surface brightness profiles for the BCG in Abell 2261. The panels (c) and (d) show the 
same for the BCG in Abell 2537. The dashed line shows the FWHM of the moffat profile seeing disk. The dotted line shows the radial 
limit to which the data are reliable. Beyond this limit the data are systematically biased by contamination from mask effects. Identical 
tests are performed on the B and R profiles for the clusters observed with CFH12k. 
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